Non-Gaussian fluctuations of mesoscopic persistent currents 
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The persistent current in an ensemble of normal-metal rings shows Gaussian distributed sample- 
to-sample fluctuations with non-Gaussian corrections, which are precursors of the transition into the 
Anderson localized regime. We here report a calculation of the leading non-Gaussian correction to 
the current autocorrelation function, which is of third order in the current. Although the third-order 
correlation function is small, inversely proportional to the dimensionless conductance g of the ring, 
the mere fact that it is nonzero is remarkable, since it is an odd moment of the current distribution. 



Already in the early days of quantum mechanics, it was 
realized that the ground state of small conducting ring- 
hke structures penetrated by a magnetic field could con- 
tain a dissipationless current circling through the struc- 
ture Such persistent currents have been well known 
in the context of superconducting devices since the early 
1960s P]. Normal-state persistent currents became the 
subject of intensive research only after it was understood 
that elastic impurity scattering, unavoidable in almost all 
realizations, reduces the phenomenon, but does not com- 
pletely suppress it [3]: The typical magnitude of a per- 
sistent current in a disordered normal-metal ring with 
circumference L and diffusion constant D is I ^ e/r^, 
where tl = L^/D is the time it takes an electron to dif- 
fuse once around the ring |4j. 

Because of time-reversal symmetry, a persistent cur- 
rent must vanish at zero field, whereas gauge invariance 
imposes that the current is periodic in the flux pene- 
trating the ring, with period 0o = h/e _5:. The precise 
value of the current I{(t>) depends strongly on the spe- 
cific disorder realization. When averaging over disorder, 
the odd Fourier components of the average current {I{<j>)) 
become vanishingly small. As long as the applied mag- 
netic field is weak (such that time-reversal symmetry is 
preserved inside the metal ring), the even Fourier com- 
ponents of {I{4>)) are nonzero, but still small compared 
to those of the typical current (/^ ((/))) |SHH]- Whereas 
early experiments reported a persistent current a few or- 
ders of magnitude larger than the theoretical prediction 
[51[TU], improved SQUID detection of the magnetic re- 
sponse of a small number of gold rings [TT] and, very 
recently, current detection through the frequency shift 
of cantilevers covered with aluminum rings |12j unam- 
biguously confirmed the predictions for the average and 
typical current amplitudes. 

The sample-to-sample distribution of I{(t>) is well ap- 
proximated by a Gaussian if the dimensionless conduc- 
tance g of the rings is large, but deviations from Gaussian 
statistics are expected to become relevant if the rings are 
very thin or strongly disordered, in which case g is small. 
The unprecedented accuracy of the cantilever method 
opens the possibility to measure the full probability dis- 
tribution of the current, not just its first two moments. 
Non-Gaussian fluctuations play an important role in the- 



oretical and experimental studies of transport in disor- 
dered metals [131 HI] , being a precursor of the transition 
to the Anderson localized regime. Non-Gaussian fluc- 
tuations were also considered in the conductance distri- 
bution of chaotic quantum dots, for which they are an 
indicator of fully quantum coherent transport |15| . 

In this Letter, we calculate the leading non-Gaussian 
connected autocorrelation function 

if('^l,02,<^3) = (/(0lW2)/(</'3))c (1) 

of persistent currents to leading order in 1/g. Here the 
brackets (. . .) denote a disorder average and the sub- 
script 'c' refers to the 'connected average', {abc)c = 
(abc) - {ab){c) - {ac){b) - (6c) (a) + 2(a)(6)(c). At equal 
values of the arguments, K gives the third cumulant 
(/((/))'^)c of the current distribution, which is propor- 
tional to its skewness. As we will show below, (/(0)^)c is 
nonzero and of order ~ e^/gr^ for small magnetic flelds 
only, and becomes vanishingly small for large magnetic 
fields. However, the full third-order correlation function 
K((j)i, (j)2, 4>3) remains nonzero for arbitrary magnetic 
field strengths. This is remarkable because K{(j)i, (j)2, (ps) 
addresses an odd moment of the current distribution. 
Previously published order-of-magnitude estimates for 
K [ini E] are a factor ~ tl /t larger than the result of 
our calculation, r being the elastic mean free time [T^. 

We will first describe the calculation of the connected 
autocorrelation function K((j)i, (j)2, (f>3) for zero temper- 
ature, at weak magnetic fields, and in the absence of 
spin-orbit scattering. Then we evaluate the effect of fi- 
nite temperature, larger magnetic fields, and significant 
spin-orbit scattering. Finally, we also discuss the possi- 
bility to identify the non-Gaussian current fluctuations 
in persistent current measurements. 

Starting point of a calculation of the persistent current 
is the thermodynamic relation [5J [TH] 

/ ^ 0) + ^{Nifi, - {N)f'^ , (2) 

which expresses the current / in a ring with a fixed num- 
ber of electrons in terms of the grand canonical potential 
and the grand canonical fluctuations of the particle- 
number N . Both J7 and N are calculated as integrals of 
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FIG. I: (Color online.) Diagrammatic representation of 
(a) the third-order autocorrelation function of the density of 
states, {i'{Ei,cf>i)i'{E2,4>2)i'{E3,cj)-i))c and of (b) the Hikami 
box T-L. The advanced (retarded) Green functions are rep- 
resented by blue dashed (red solid) lines. The thick gray 
lines connecting the Green functions represent the diffuson 
and cooperon ladders. The dotted black lines represent im- 
purity scattering. 



autocorrelation function (i/(E'i, 0i)z^(i?2, </'2)'^(-E'3, 03))c. 
Performing the disorder average using diagrammatic per- 
turbation theory [20], the latter can be represented 
schematically as in Fig. [T^. In this figure, the advanced 
(blue dashed lines) and retarded (red solid lines) Green 
functions are connected to each other by thick gray lines 
which represent the diffuson (— ) and cooperon (-f ) lad- 
ders, (Qq^; Eap), where we use the short-hand no- 
tations q^^ — Qq, ± and Eo,p — — Efj. In this 
notation, the index a refers to the wave vector, energy 
and spin of the advanced Green function and (3 to those of 
the retarded function. In the absence of a magnetic field 
penetrating the metal ring and without spin-orbit scat- 
tering, the diffuson and cooperon propagators are given 
by the standard expression 



^{vWD\ci%\^ + ^E^0- 



(5) 



the density of states v{E, 
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where T is the temperature and the chemical potential, 
which is chosen such that the disorder average (iV(/i)) 
equals the (canonical) particle number N . 

Since r2(/i, (f) and A^(/x, (j)) are statistically uncorrelated 
to leading order in l/g — they have opposite parity un- 
der a particle-hole conjugation — , only the first term in 
Eq. ([2]) contributes to the autocorrelation function K. 
Calculation of K then proceeds through the standard 
relation between the density of states J^(£', (j)) and the 
Green function, 
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where a is the spin index. 

Calculation of the current autocorrelation function K 
requires calculation of the connected density-of-states 



The ladders are connected at a 'Hikami box' % depicted 
in Fig. [1)3, 
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where the indices s and s' can be -I- or — , depending 
on whether the upper (s) or lower (s') ladders repre- 
sent diffuson or cooperon propagators. Combining the 
expressions for T) and "H, one then computes the third- 
order density-of-states autocorrelation function by sum- 
ming over the momenta qi , q2 , and qa , the corresponding 
spin indices, and the diffuson/cooperon indices s and s' . 

We now take the ring's circumference L to be much 
larger than its thickness. In that limit, only terms in 
which the transverse component of the momenta q* van- 
ishes contribute to the correlation function. In the pres- 
ence of a fiux (j)a penetrating the ring, the longitudinal 
component q of the momentum takes the discrete values 



qa = - '■Pa), With 
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n being an integer. Converting the sums over the integers 
n to integrals using the Poisson summation rule, we then 
find that, at temperature T = 0, 



^(^1' <^2' ^3) = ^ ^ ^ ^ j:^^ j:^ j:^ M^^k^.,) M^^k ^^,), (s) 

I 

where the dimensionless conductance g = 2Trh{iy) /tl, and with the notation ip^g = ipazL Lpp. The scale of the 
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non-Gaussian fluctuations, when normalized to the Gaus- 
sian fluctuations, is K^/'^ / {l'^)^/'^ ~ g~^^^- The magni- 
tude of the summand in Eq. (|8| drops to zero quickly for 
increasing k and k' , so that only the lowest harmonics 
will contribute signiflcantly to K. Note that the third 
cumulant. 
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fc2 + 4fcfc' + 

k/{k + k')^ (9) 

X sin(47rfciy9) cos(47rfc'9?), 



is non-zero only due to the 'double cooperon' contribu- 
tion with s ^ s' ^ +. 

Let us now discuss the effect of a finite magnetic field, 
spin, and finite temperature. All three effects are relevant 
for a realistic description of the experiment of Ref. [H] . 
Following Ref. [21] , we model the magnetic field penetrat- 
ing the sample as toroidal. Although a toroidal magnetic 
field does not fully correspond to the experimental ge- 
ometry, it is believed to provide a reliable and tractable 
description of the experiment ^21J . With a toroidal mag- 
netic field -By , the transverse component of the momenta 
is no longer zero, which leads to a shift of the mo- 
mentum eigenvalues appearing in the expression for the 
cooperon ladders. 



(10) 



e_L being a geometry-dependent positive constant propor- 
tional to . For a ring with circular cross section with 
radius R, one has e\ = {l/^TT'^h'^){eLRB\\f + 0[Bf^ [H]- 
When the coupling of the applied magnetic field to 
the electron spin (with Zeeman frequency uji) or spin- 
orbit scattering (with spin-orbit scattering time Tgo) are 
included, the diffuson and cooperon propagators acquire 
additional structure with respect to their spin indices. 
As long as t/tso ^ 1, the Hikami box remains spin- 
conserving, such that elements of the propagators which 



flip the spin of one of the Green functions do not con- 
tribute. In this case only 'diagonal' elements play a role 
and the two spin indices in Eq. JsJ continue to be suf- 
ficient. With the toroidal field (10 1 included, the spin- 
dependent propagators read 
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where I? (q , E) is the spin- independent propagator of 
Eq. (Isl, and the factors a;^ ^ ^ ^""^ defined as 
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For a large Zeeman splitting {tlujz 3> 1), the contri- 
bution of the elements 2^^,^ and 2?^,^ will be suppressed, 
which leads to a factor four decrease of K{(j)i,(j)2, (^3) with 
respect to the zero-field case of Eq. ([8| . Strong spin-orbit 
scattering {tl/tso ^ 1) causes an additional factor four 
decrease of K{(j)i,(j)2, ^3). 

Inserting the diffuson and cooperon propagators of Eq. 



( 11 1 and extending the previous calculation to finite tem- 



peratures, we arrive at the complete result 
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(13) 



where 9 — 2T:kBTTL/h is the rescaled temperature. 



Equation ( 13 ) shows explicitly how the various contributions to K are suppressed by a magnetic field, spin-orbit 
scattering, and finite temperature. Using characteristic values for the experiment of Ref. |12j . L ~ 2.5 fim, R ^ 50 nm 
and ~ 3.5 T, one finds that e± ~ 12. Since all x'^ ^ contain a term An'^e^, a large magnetic field, e\ ^ 1, 
annihilates the contribution of all cooperons, and thereby also the third cumulant of Eq. In this regime, one 
thus indeed expects to find (/(0)'^)c ~ 0. However, since x^iq_i = irrespective of the magnetic field and spin- 
orbit scattering strengths, the general connected correlation function remains nonzero at high fields and with strong 
spin-orbit scattering because of the contribution from the diffusons and 'D^ i^- 

A large temperature leads to exponential suppression of all contributions to the connected correlation function 
K((j)i, (j)2, (1)2) and to a suppression of the non-Gaussian fluctuations in comparison to the Gaussian fluctuations. In 
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the high-temperature Umit, 0^1, and for small magnetic fields and weak spin-orbit scattering, tl/tsotTlujz ^ 1, 
(but still assuming e\ ^ 1), Eq. (13 1 simplifies to 



I sm[2n{ipi +ip2- 2(^3)] + sin[27r(^2 + ^3- 2(pi)] + sin[27r(933 + (fi - 2(^2)]}- (14) 



For comparison, in the same limit the second-order cur- 
rent correlator, which describes the Gaussian fluctua- 
tions, scales as {P) oc 9^e~"^ jB], so that typically 
^i/3/^j2^i/2 _ g-v^/6^-1/3 investigate the 

full temperature dependence of K, we show in Fig. [2] the 
temperature dependence of the leading Fourier compo- 



nent Ki,i (which is the one given in Eq. (14) above) 



As a complete measurement of the connected correla- 
tion function A'(0i, 02, '/'s) is likely to be very involved, 
it is worthwhile to investigate the behavior of the cor- 
relation function when two of its arguments are equal, 
K((j), (f), (f) + A(f)). In the same limit of small magnetic 
fields and weak spin-orbit scattering, but with 3> 1, 
this correlator is a function only of the flux difference, 
K{(p, (pjcf) + A(j)) — K'{Aip). The zero temperature limit 
for K'{A(/)) is shown in the inset of Fig. [2] In the experi- 
ment of Ref. [12] one has roughly g ~ 10 for the smallest 
rings, so that ^K'{A(j)) ~ 0.1 nA. We thus believe that 
a measurement of K'{A(p) should be feasible. 

In conclusion, we have evaluated the third-order au- 
tocorrelation function for the persistent current in an 
ensemble of disordered metal rings. Remarkably, this 
correlator remains non-zero at large magnetic fields, al- 
though the non-Gaussian current fiuctuations are a factor 
~ g^/^ smaller than the leading Gaussian fluctuations. 
The third-order connected correlator is the leading-order 
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FIG. 2; (Color online.) (Main) The Fourier component 7^1,1 
as a function of the rescaled temperature 8 — 2nkBTTi^/h,. 
Ki,i is normalized to the zero-temperature value predicted by 
Eq. pi For > 25 the high-temperature approximation given 
by Eq. 14 (red solid line) and a numerical evaluation using Eq. 
13 (circles) agree well. (Inset) The correlator K{(f), 4>, -I- A^) 
at zero temperature in units of e^/vr^rff;. For this plot we 
used Eq. [s] summing over fc, k' < 40. 



signature of localization effects. About one decade ago, 
the observation of non-Gaussian correlations in the trans- 
port of disordered quantum wires and chaotic dots was 
greeted as a unequivocal signature of quantum coherent 
transport [I3HT5| . We hope that our calculations moti- 
vate for a corresponding breakthrough in the realm of 
persistent currents and other equilibrium properties. 
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